Theory of Nonlinear Dispersive Waves and Selection of the Ground State 
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A theory of time dependent nonlinear dispersive equations of the Schrodinger / Gross-Pitaevskii 
and Hartree type is developed. The short, intermediate and large time behavior is found, by de- 
riving nonlinear Master equations (NLME), governing the evolution of the mode powers, and by a 
novel multi-time scale analysis of these equations. The scattering theory is developed and coherent 
resonance phenomena and associated lifetimes are derived. Applications include BEC large time 
dynamics and nonlinear optical systems. The theory reveals a nonlinear transition phenomenon, 
"selection of the ground state", and NLME predicts the decay of excited state, with half its en- 
ergy transferred to the ground state and half to radiation modes. Our results predict the recent 
experimental observations of Mandelik et. al. 11} in nonlinear optical waveguides. 



Recent exp erimental [H, 0, 0, and theoretical 
|1 US 003 breakthroughs in the field of Bose- Einstein 
Condensates (BEC) and the intense interest in nonlin- 
ear optical devices 01 necessitates the detailed analysis 
of time-dependent nonlinear dispersive equations, such 
as the time-dependent Gross-Pitaevski (GP) and nonlin- 
ear Schrodinger (NLS) equations in one or more spatial 
dimensions. Also, new coherent phenomena in soft con- 
densed matter, such as macroscopic resonances 0, 124^ . 
DNA-denaturation dynamics and macromolecule dynam- 
ics, have rekindled the need to determine the large time 
behavior of time-dependent Hartree - Fock (HF) type 
equations. In contrast to linear scattering, where time- 
independent methods are used to compute approximate 
S'-matrix for general multichannel processes, nonlinear 
problems are not amenable to such theory, except in com- 
pletely integrable cases. Thus, mainly large scale numer- 
ical methods have been applied so far. 

In this Letter we describe and apply a rigorous time- 
dependent theory and approach developed in to pro- 
vide detailed explicit behavior for a class of nonlinear 
dispersive partial differential equations, which include 
NLS, GP and HF, on short, intermediate and infinite 
time scales. We show, as a consequence of our analysis, 
that in a multimode nonlinear dispersive systems which 
has multiple nonlinear bound states, that nonlinear reso- 
nant interactions lead to a crystallization of the coherent 
and localized part of the solution on the nonlinear ground 
state (selection of the ground state); the nonlinear excited 
states are metastable and decay at a rate given by a non- 
linear analogue of Fermi's golden rule. 

In very recent experimental work pdj on CW beams 
centered at wavelength, A = 2ir/ko, propagating in mul- 
timode nonlinear AlGaAs waveguides, the output distri- 
bution of optical power is measured, and ground state 
selection and the partition law Q} have been demon- 
strated. Our results predict these observations and can 
be used to interpret the earlier work on photorefractive 
waveguides [T^ . 

For definiteness we consider the nonlinear Schrodinger 



equation in three space dimensions: 

id t <j> = H cj> + gF(\<f>\ 2 )<p. 
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F denotes the 

m\ 2 ) = h 2 , 

type type F(\tp\ 



nonlincarity, here taken to be local, 
but more generally, of Hartree - Fock 
! ) = J G(x,y)\ip(y)\ 2 dy with, for exam- 
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> 0. g denotes the 



pie, G(x,y) = \x 
coupling coefficient, related to the scattering length in 
BEC and, in optics, to — 712, the nonlinear Kerr coeffi- 
cient. Hq = —A + V(x), where V(x) denotes a poten- 
tial which decays to zero sufficiently rapidly at infinity, 
and for which Hq has two bound states (for simplicity); 
Hoipj* = Ej*ipj*, j — 0,1. Such equations appear in 
many applications in which coherent soliton-like struc- 
tures and their internal mod es l22| interact with disper- 
sive waves, e.g. [Jlll!E0|2iTThe effective equations 
of an atom in a field or active media is also of this form 
0,113 ■ The requirement on V(x) that H have (at least) 
two bound states makes the model of general interest and 
the phenomena very rich. 

Nonlinear bound state solutions bifurcate from linear 
bound states for weakly nonlinear perturbations |14| giv- 
ing rise to nonlinear ground state family and nonlinear 
excited state families (for two bound state Hamiltoni- 
ans). Specifically, for energies close to i?o* and J5i*, 
the eigenvalues of Hq, there are nonlinear bound states: 
e~ iEot ipE , ^E x er lExt , which solve QJ. i\>E a and ipEx 
are complex-valued (exponentially) localized solutions of 
(H + gF(\4, E% \ 2 )) $ Ei = EiipEi ■ When g = 0, E = £ * 
and Ei — Ei* are values for which there are nontrivial 
solutions. For j ^ there is a solution for all E near 
Ei*. ipEi can be thought of as a pinned soliton, ground 
(J = 0) or excited state (j = 1), within the support of 
V(x). We assume w* = 2Ei* — E * > 0, ensuring cou- 
pling of bound states to the continuum at second order; 
analogous results with slower decay rates hold more gen- 
erally. 

The main results we establish concern the time evolu- 
tion of solutions to QJ for all data (with small energy). 
The large time behavior is as follows: For any initial data, 
with small energy the solution ip(t) of QJ is asymptoti- 
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cally given by 

m - 
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e ip±, 



as t — > ±00, 



in L 2 . Generically, j = 0, corresponding to the ground 
state, while there is a finite dimensional set of data for 
which j = 1. Thus, for typical initial conditions we have 
ground state selection. We analyze the detailed dynamics 
by reduction of the NLS to a finite dimensional system, 
which dominates the full dynamics. 

The main second result concerns the transient, finite 
time, behavior of the system, which is very rich: There 
are three time scales, determined by the initial data. In 
the first time interval, the dynamics is dominated by ra- 
diative terms (formative stage); once the initial radiative 
part propagates far enough from the support of V(x), the 
soliton's position, the second period (embryonic stage) 
begins. This is marked by a monotonic (exponentially 
fast) increase in the ground state amplitude, relative to 
the excited state amplitude. This stage proceeds as long 
as the excited state part is larger than some fraction of 
the ground state part. Once the ground state reaches 
a size comparable with some (fixed) fraction of the ex- 
cited state, the third and final stage begins (ground state 
selection); the ground state amplitude increases mono- 
tonically and the excited state decreases monotonically, 
but at a polynomial rate as t — > 00. (This describes 
the generic case; the asymptotic state is given by an ex- 
cited state a 2(| whenever the first stage persists to time 
equal infinity.) A different approach, based on lineariza- 
tion around the excited state for intermediate times and 
around the ground state for large times, gives similar re- 
sults (Soffer-Weinstein, unpublished notes, 1996). 

This phenomena is quantified by our derivation of the 
Nonlinear Master equations (NLME) that govern the dy- 
namics of the coupled ground and excited state modes 
(and the radiation modes . . .). If we denote by Po(t) and 
Pi (t) , respectively, the (up to near identity transforma- 
tions) squared projection of the system's state onto the 
ground state and excited states, at time t, we have 



d t P = 2rP P 2 + pa(t)PaP? I( > + 0(r 5 ' 2 ) 

dti\ = -4i\p p 2 + Pl (t)y/p~ p^ + e>(£- 5/2 ), 

r = ng 2 I T [dF[ipo*i/ju]ilJu] (w*) | 2 
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dF[lp *lpl*}^l 



G(x - y)ipo*(y)ipu(y)dy ipi*(x). 

(.3) 



PM (t) and pMi (t) are small oscillatory functions of time, 
Mo, Mi > 3 and T[q](oj) denotes the projection of q 
onto the generalized eigenfunction of Ho at frequency 
lu. The crucial number T is (generically) positive and 
given explicitly in terms of known eigenstates of Hq. 
r ^ is the nonlinear Fermi Golden rule for such sys- 
tems and it gives the rate of decoherence and relaxation 
[Hill El El i3 . The behavior of these Master Equa- 
tions (£2J) reflects the three time scales mentioned above, 



on which the behavior is very different. For large enough 
time, in the third time domain, the last two terms in 
can be ignored. In this latter regime, it is easy to see that 
Pi(t) -> as t -> 00 and d t (2P + Pi) = 0. Therefore, 
2P (oo) = 2P (t ) + Pi (to)- K follows that 



P (oo) = P (t ) + \Pi(t ); 



(•4) 



half of the excited state energy flows to the ground 
state (the other half goes to radiation), a kind of en- 
ergy equipartition. (The factor 2 also appears in the 
linear analysis, and is interpreted as the ratio of re- 
laxation to decoherence time [2! |25|.) Furthermore, 
it follows that for any initial state (Po(fo), Pi(io)) with 
Po(to) 7^ the system converges to (Po(oo),0) with a 
rate ~ (1 + 4I\Po(oo)Pi (to)t)^ 1 (selection of the ground 
state). Combining the analysis of the NLME Jj2J for 
all time scales with the previous statements gives a com- 
plete description of the solution for all time scales and, in 
particular, the asymptotic behavior, relaxation and de- 
coherence rates, the asymptotic profile and energy of the 
soliton/ground state. 

We sketch our method, for the special case of QJ 
with cubic nonlinearity, F(\(j)\ )<j> = \< t»\ 2 (( >- Our approach 
makes use of ideas from pjl ^1 IT7L fl8L fis| . We begin 
with the Ansatz <j>{t) = e~ ie ^[tp (t)+tpi(t)+(t)2(t)} where 
i>o(t) = 4>E a (t) is a solution of the ground state eigenvalue 
equation with energy Eo(t), at time t. Eo(t) will be de- 
termined later by orthogonality conditions 0, 0, . 
Similarly ipi (t) is an excited state eigenvector with eigen- 
value Ei(t). 9(f) = e (t)+8(t), 6 (t) = J Q t E (s)ds. 
6(f) will be chosen appropriately; it includes (logarith- 
mic) divergent phase. Substitution of the above Ansatz 
for (p into ifTTl) . and complexifying the equations (4>2 — * 
(02, <h) = $a(*)) (V>j (*Pj,~4>j) = et c-) we derive 

id t <5> 2 (t) =Ho(t)<S> 2 (t) - idt^o 

((E Q - Ei) + d t 8)a 3 + id t ] *i+F NL , (.5) 

where P/vl is nonlinear in $2, ^0, ^1, 8 and TCo(t) is 
given by the matrix operator 



0-3 



' H - E (t) + 2g\Mt)\ 



9*Po(t) 
H-E (t) + 2g\Mt)\ 



where 03 is the Pauli matrix diag(l, — 1). We consider 
the spectrum of TLo(t) for fixed t and l^ol = l a o| small: 
(a) The continuous spectrum extends from — fi to —00, 
and \i to 00 where \i = Ei — Eq + O(|ao| 2 ). The dis- 
crete spectrum is {0,— with < < \Eq\ by 
assumption, (b) Zero is a generalized eigenvalue of Ho, 
with generalized eigenspace. spanned by {a^^Oi dE a ^o}- 
The discrete spectral subspace has dimension four. 
Therefore, $2 which lies in the continuous spectral part 
of Ho(t), is constrained by four orthogonality condi- 
tions. Furthermore, dtO is chosen to remove divergent 
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logarithmic phase contributions. In the weakly non- 
linear (perturbative) regime, bound states have expan- 
sions ipEj = aj(ipj*(x) + g\aj\ 2 i/)f\x) + 0(g 2 \aj\ i )) and 
Eg = Ej*+0(\aj\ 2 ). The system for $ 2 and a = (do, «i) 
can be written in the form: idtd — A(t)a+F a , idt&2 = 
H(t)<£ 2 + F<i> . 

To study the energy exchange between ao and a%, 
it is important to express ai(t) as a slow amplitude 
modulation of a rapidly varying phase. With this goal 
in mind, we study the equation for ai, expressible as 
idtdi = A(t)di + F. Freezing t at some arbitrary large 
time T in A(t), we solve idtdi = A(T)a\ which gives 
a periodic Floquet solution matrix X(i). We use that 
to eliminate this (fast) oscillation in d\ = X(t)(3%. The 
difference A(t) — A(T) and other similar differences are 
higher order corrections, uniformly in t, as t, T — » oo, 
and therefore can be neglected. 

To proceed further we decompose $2 into its continu- 
ous spectral (dispersive) part, 7/, relative to Ho(T), and 
its components along the discrete modes. The latter are 
higher order and controllable. Thus NLS, at low energy 
is equivalent to a system of the form: 

id t r) = H (T)n + F n (t; a , 
idA = 2g(^^l)\f] 1 \ 2 aoe^- E ^ t 

+ giMl,^)^- 2 ^-^ + Ki, 

where IZj denotes corrections of a similar form and higher 
order. 

The above system can be viewed as an infinite dimen- 
sional Hamiltonian system consisting of two subsystems: 
a finite dimensional subsystem governing "oscillators" , 
(ao , fli ) , and an infinite dimensional subsystem governing 
the field, r\. Although this system has time-rapidly vary- 
ing coefficients and no evident direction of energy flow, we 
claim we have now made explicit the key aspects, which 
give rise to resonant energy exchange. This is made ex- 
plicit via a detailed analysis of non-resonant and resonant 
terms Non-resonant oscillatory dependence contributions 
can be transformed by successive near-identity changes 
of variables, (ao,/?i) <— > (<So,/3i), to higher order order 
in the energy of the data (assumed small) and perturba- 
tively controlled. Resonant oscillatory terms cannot be 
transformed to higher order and contribute to the finite 
dimensional reduction. 

To arrive at the reduction, we solve the r\ equation, 
making explicit all terms through second order in g, using 
the Green's function G(t,t') = e -™-o{T)(t-t') _ We focus 
on the key terms coming from the sources in T n or the 
type a a\, < i,j < 2 and having oscillatory phases 
e*""'*. Their contribution to 77 is of the form 

~ / e- m °( T )^- t ')|x>e im « (i ' ) 4(t / )ai(t / )^ / (-6) 
Jo 



where ao,ai is a component of either do or c?i, where 
x) is an (exponentially localized) function of position, 
expressible in terms of ipo* and ipi*. We insert this solu- 
tion into the ao, a\ equations, in place of $2- We obtain 
integro-differential equations for ao, ai,(/3i). Terms of 
the form are solutions to a forced linear system and 
among the forcing terms made explicit in (£5J) are oscil- 
latory terms with the frequency w*, which are resonant 
with the continuous spectrum. Internal dissipation re- 
sulting in nonlinear resonant energy transfer from the 
excited state to the ground state and to dispersive ra- 
diation is derived from these resonant terms; see also 
the derivation of internal dissipation in both linear and 
nonlinear resonance theories, recently developed by us 
El 111 El E3 This dissipation coefficient is T, the 
rate of decoherence and relaxation. The above described 
scheme gives id t do = (— A + iT)\Pi\ 4 do + TZo(t), id t $i = 
2(A~iT)\d \ 2 \Pi\ 2 (3i+n 1 (t). 

Introducing the squared projections of the system's 
state onto the ground state and excited states, Po = 
\d \ 2 , Pi = \Pi\ 2 we obtain NLME, ©■ The system 
(EU is analyzed in terms of renormalized powers, Qo 
and Q\, for which it is shown that there exist transi- 
tion times to and t\, such that: Qo(t) decays rapidly on 
[0, to], Qo(t)/Qi(t) grows rapidly on [toj^i] and the fi- 
nally on [ti,oo) the following system governs: dtQo = 
2TQoQ\, d t Qi = -4TQ Qi- This gives Q T Qo(oo) 
and Qi | at rates discussed above. 

The decay of the nonlinear excited state can also be 
understood as a linear instability. Let TL\ denote the lin- 
earization about the nonlinear excited state, the operator 
in (J3J with tpo replaced by the excited state ipi. Since 
cj* > 0, in the limit of vanishing nonlinear terms, Hi 
has an embedded discrete eigenvalue in the continuous 
spectrum. Perturbation theory of embedded eigenvalues 
(see, for example, the time-dependent approach in [l5| 
), can be applied to show that this embedded eigenvalue 
perturbs to a linear exponential instability with expo- 
nential rate ~ T associated with the linear propagator 
exp(— iHit). 

In the context of nonlinear optical waveguide experi- 
ments, e.g. our analysis gives the distance, Ltransf, 
over which nonlinear resonant energy transfer occurs as: 
L t~ra„sf = 4rPo(oo)P x (0). Here g = -n 2 /n , n 2 [m 2 /W] 
the Kerr nonlinear coefficient, Pj [V 2 /m 2 ] the square of 
the electric field projection onto state j. T is obtained 
from Q5J), where the states ?pj* are those derived from 
the Hamiltonian, Ho = — 2kg 2 A x + V(x), with potential 
V(x) — (1/2)(1 — (no(x)/fio) 2 ), with x, the transverse 
coordinate. no(x) and no are, respectively, the linear 
refractive index and its value at infinity. This implies 

Klsf - (4*** (Jt*) 2 • fe) 2 • |^w&]K)| a ). 

For parameter values at wavelengths ~ 1.5/uto and peak 
power levels at I0 3 W, i trans f of device size dimensions 
(100 fim to 1 mm) can be attained. The factor ^[^(o;*), 
which can be approximated by WKB, depends on the 
density of states of Ho near a;*, and can be tuned by 
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varying the waveguide material and geometric parame- 
ters. A second application to optical devices is the use 
of periodic photonic microstructures with appropriately 
designed defects to trap coherent light pulses. For ex- 
ample, it has been shown theoretically that (Kerr) non- 
linear and periodic structures support gap solitons pflj 
and trap pulses [3(j traveling at any speed less than the 
speed of light. Experiments have demonstrated soliton 
propagation at about 50%c and theoretical studies 
[2l| show that gap solitons can be trapped with appropri- 
ately designed defects. This has potential applications to 
optical buffering, high-density storage and optical gates. 
The light stopping mechanism, is based on the transfer 
of energy from incoming solitons to the pinned nonlinear 
modes of the defect. The ground state selection phe- 
nomenon and the partition of excited state energy into 
pinned nonlinear ground state and radiation modes QJ), 
described by NLME, quantify the rate of transfer and 
efficiency of trapping. For systems of BEC droplets in 
a potential well, the relaxation time to the ground state 
and the decoherence time due to the presence external 
perturbations, e.g. other droplets, are measurable and 
of fundamental importance related to T. For example, 
they are relevant to the construction and feasibility of 
quantum gates and memories. 

Our analysis is applicable in Hamiltonian systems 
which can be decomposed into a lower dimensional dy- 
namical system (oscillators) coupled to an infinite dimen- 



sional dynamical system (wave- field), acting as a disper- 
sive "heat bath", e.g. how an open quantum system is 
effected by its environment. T is related to the rate of 
energy transfer or decoherence. In general, when a mean- 
field type approximation to the interaction with the en- 
vironment describes the system, the above analysis can 
be used, e.g. to study wave function collapse in various 
systems. 

To summarize, we have derived the behavior on all 
short, intermediate and infinite time scales of NLS type 
equations, in which general solutions involve an inter- 
action among two families of nonlinear bound states 
(soliton-type) and dispersive radiation. We have derived 
Nonlinear Master equations, which govern the dynamics. 
The phenomena of ground state selection was demon- 
strated, the rate of decay /relaxation was computed, and 
the energy distribution of the asymptotic state of the 
system was derived. These phenomena were recently ob- 
served experimentally in [11]. Finally, we have shown 
the applicability of our theoretical results to nonlinear 
optical devices and to BEC dynamics and decoherence 
phenomena. 
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